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1.  Introduction 

In  recent  years,  there  has  been  a  great  deal  of  interest  in  shock  propa¬ 
gation  in  non-uniform  ducts.  Chester  [1]  determined  the  disturbance  produced 
behind  a  shock  of  arbitrary  strength  propagating  through  a  tube  of  arbitrarily 
varying  cross-section  by  linearizing  the  problem  on  the  basis  of  small  area 
variations  and  found  that  the  pressure  perturbation  behind  the  shock  was  given 
by 

-  K(P2  -  Pi)[AA]/A  , 

where  P  -  P.  was  the  initial  pressure  discontinuity  across  the  shock,  [AA] 

w  X 

the  net  change  in  area  and  the  parameter  K  a  monotonically  decreasing  function 
of  the  shock  strength.  Chester  started  with  the  full  three-dimensional  equations 
of  motion  and  then  carried  out  an  averaging  process  by  restricting  the  final 
consideration  to  the  average  pressure.  However,  such  a  restriction  at  the  end 
means  that  the  same  results  must  be  obtainable  by  starting  with  a  one-dimensional 
analysis.  This  important  fact  was  first  realized  by  Paul  Germain,  and  Chester's 
results  obtained  by  a  one-dimensional  analysis  by  Gundersen  [6],  extracts 
from  which  appear  in  [7],  the  analysis  being  based  on  techniques  presented  by 

*  Sponsored  by  the  Mathematics  Research  Center,  United  States  Army,  Madison, 
Wisconsin  under  Contract  No.  DA-ll-0<iZ-ORD-2059. 
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Germain  and  Gundersen  [5].  Surprisingly,  the  one-dimensional  approach  leads 
to  an  improvement,  namely,  the  term  [AA]/A  in  Chester' s  work  is  replaced 
by  A/A  where  A  is  the  perturbed  area  distribution.  This  shows  the  relation 
between  the  area  change  and  shock  strength  at  any  point.  This  fact  was 
observed  by  Chisnell  [2], who  obtained  it  in  a  different  manner,  essentially 
assuming  Chester's  steady  state  solution, valid  for  large  time  only,  could  be 
utilized,  and  used  to  apply  an  integrated  form  of  Chester's  solution  to  discuss 
converging  cylindrical  and  spherical  shocks  with  great  accuracy. 

Lately,  many  others  have  worked  on  similar  or  related  problems,  e.g., 
Rosciszewewski  [8]  and  Whitham  [9]  who  also  rederived  Chester's  results. 

In  the  present  paper,  the  techniques  presented  in  [7]  are  generalized  to 
determine  the  perturbation  produced  when  an  initially  plane  magnetohydrodynamic 
shock  wave  encounters  an  area  variation,  the  problem  being  linearized  on  the 
basis  of  small  area  variations.  In  order  to  ensure  that  all  changes  in  the 
behavior  of  the  shock  are  caused  by  the  non-uniform  area,  it  is  assumed  that 
the  area  variations  are  confined  to  the  region  x  >  0  while  to  the  left  of  the 
section  x  =  0  (say),  the  tube  is  of  constant  area,  and  the  initially  plane 
shock  propagates  through  this  portion  with  constant  speed.  The  fluid  in  front 
of  the  shock  is  assumed  at  rest. 

When  the  shock  meets  the  area  variation,  it  is  perturbed,  the  shock  strength 
altered,  and  the  subsequent  flow  is  non-isentropic. 

There  are  two  distinct  contributions  to  the  disturbance,  namely,  a  per¬ 
manent  perturbation  caused  directly  by  the  area  change  and  a  transient 
disturbance,  due  to  reflection  from  the  shock  of  the  permanent  perturbation. 
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2  2  1/2 

which  propagates  with  velocity  w  =  [c  +  b  ]  ,  where  c  is  the  local  speed 

of  sound  and  b  the  Alfven  speed,  with  respect  to  the  flow  behind  the  shock. 

In  the  immediate  vicinity  of  the  shock,  the  ultimate  effect  is  an  altered 
shock  strength  and  concomitant  pressure  change  behind  the  shock.  When  the 
flow  behind  the  shock  is  <  a>,  there  is  a  disturbance  convected  to  the  left  with 
velocity  q  relative  to  the  fluid  while  for  flow  >  u,  a  similar  disturbance  is 
generated  which  propagates  to  the  right  with  velocity  u  relative  to  the  fluid. 
This  latter  must  be  added  to  the  steady  flow  solution.  Expressions  for  these 
various  contributions  are  obtained. 

Specifically,  it  is  found  that  the  pressure  perturbation  immediately  behind 
the  incident  shock  is  given  by: 

=  ~  mj^)(P2  - 

where  perturbations  are  denoted  by  a  bar,  P^  -  P^  is  the  original  pressure 
discontinuity  across  the  shock,  the  cross-sectional  area,  tr  the  shock 
strength  defined  as  a  density  ratio  and  m^^  =  bj^/c^^  a  measure  of  the  applied 
field. 

For  all  mp  lim  K-  =  0.5  and  lim  =  0.45085  [the  present  analysis 
0—1+  0--4- 

is  limited  to  a  monatomic  gas  so  that  1  <  c  <  4] .  But  these  are  exactly  the 
limits  for  the  ordinary  gas  dynamic  case  which  corresponds  to  m^^  =  0  .  Thus, 
the  result:  for  very  weak  or  very  strong  shocks  [in  an  asymptotic  sense],  the 
results  are  independent  of  the  magnetic  field  and  agree  with  the  usual  gas 
dynamic  results. 

In  magnetohydrodynamics,  strong  shocks  occur  for  o-  close  to  4  or  for 


a  very  strong  applied  field  for  any  o-  >  1 ,  i.  e. ,  is  then  large. 

For  nij  =  0,  is  a  monotonically  decreasing  function  of  the  shock 
strength  and  agrees  exactly  with  the  corresponding  parameter  in  the  afore¬ 
mentioned  publications. 

For  any  m^  ^  0,  the  monotoncity  is  lost  but  each  curve  is  concave  upward 
with  curves  for  greater  m^  lying  below  those  for  lesser  m^,  and  all  curves  pass 
through  the  points  (a,  K^)  =  (1,0.5)  and  (4,0.45085). 

Further,  for  fixed  incident  shock  strength,  is  a  monotonically 
decreasing  function  of  mj^,  i.  e.,  in  a  diverging  (converging)  channel,  the 
pressure  decrement  (increment)  is  decreased  (decreased)  by  increasing  the  applied 
field. 

If  the  area  variations  are  confined  to  a  transition  section  of  finite  extent 
joining  two  portions  of  constant  area,  the  total  effect  of  the  passage  through 
the  transition  section  is  a  change  in  shock  strength  with  the  shock  eventually 
again  becoming  uniform. 

Qualitatively,  the  motion  of  the  shock  is  independent  of  whether  the  main 
flow  behind  the  shock  is  <  w  or  >  w. 

Since  the  analysis,  which  is  limited  to  a  monatomic  gas,  includes  the 
usual  gas-dynamic  results  as  a  special  case,  there  is  a  check  on  the  theory 


presented  herein. 
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2.  General  Theory 

The  qua  si-one-dimensional  non-steady  flow  of  an  ideal,  inviscid,  per¬ 
fectly  conducting  monatomic  compressible  fluid,  subjected  to  a  transverse 
magnetic  field,  i.e.,  the  Induction  B  =  (0,0,  B),  is  governed  by  the  system  of 
equations: 


+  uc^  +  cu  j/^3  +  ucA ^3A  =  0 

(1) 

2  y  2  , 

3cc^  ^  ® ~  ^ 

(2) 

B^  +  uB  +  Bu  =0 

t  X  X 

(3) 

■H  '  0 

(4) 

7.  7 

where  u, c, t], b  =  B  /h-p  ,  p,p  and  A  are,  respectively,  the  particle  velocity, 
local  speed  of  sound,  specific  entropy,  square  of  the  Alfven  speed,  density, 
permeability  and  cross-sectional  area.  Partial  derivatives  are  denoted  by 
subscripts  and  all  dependent  variables  are  functions  of  x  and  t  alone  save 
A  which  is  considered  time  independent.  The  characteristics  of  this  system 
are: 

cbo^dt  =  u,  u,  u  +  CO,  u  -  CO 
2  2  1/2 

where  co  =  [c  +  b  ]  which  corresponds  to  the  limiting  case  of  a  fast  wave. 

For  an  arbitrary  isentropic  constant  area  flow,  the  base  flow  in  the 
neighborhood  of  which  perturbations  will  be  considered,  the  basic  system  of 
equations  (1),  (2)  and  (3)  may  be  written  in  the  following  characteristic  form 

[31: 
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Xp  -  (u  +  a))tp  =  0 

(5) 

X  -  (u  •'  w)t  =  0 
a  a 

(6) 

Xt  -  ut.  =  0 

(7) 

2  2 

(w  -  c  )Bp/B  +  uUp  +  3cCp  =  0 

(8) 

(w^  -  c^  )B  /B  -  oju  +  3cc  =  0 
a  a  a 

(9) 

B./B  -  3c. /c  =  0 
%  § 

(10) 

with  characteristic  parameters  (a,  p,  |). 

3 

From  equation  (10),  it  is  clear  that  the  quantity  B/c  is  constant  along 
each  particle  path,  which  is  a  well-known  consequence  of  the  assumption  of 
infinite  electrical  conductivity,  i.  e.,  the  magnetic  field  is  "frozen"  into  the 

3 

fluid.  For  a  cons^-ant  state  or  a  simple  wave  flow,  B/c  is  constant  through¬ 
out  the  flow,  and  equations  (8)  and  (9)  may  be  written  as: 


Up  +  3a)Cp/c  =  0 

-  u  +  3«c  /c  =  0 
a  a 


Since 


B  =  r^^c' 


p  =  r^c  with  constants  r^  and  r^  , 
=  r^^cV^r^  =  kc  ^  (say) 


so  that 


0)^  =  c^(l  +  kc) 


(11) 

(12) 


and  (11)  and  (12)  may  be  integrated  explicitly  to  yield: 
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a/2  +  (1  +  =  u/2  +  (w/c)Vk  =  a  (13) 

-  u/2  +  (1  +  =  -  u/2  +  (u/c)Vk  =  P  (14) 

where  (a,  p)  may  be  considered  as  generalizations  of  the  usual  Riemann  in¬ 
variants. 

It  is  convenient  to  make  a  transformation  of  dependent  variables  so  that 
equations  (13)  and  (14)  become  linear  relations  in  order  to  bring  the  present 
results  into  close  analogy  with  ordinary  gas-dynamics.  This  is  effected  by  the 

3 

substitution  w  =  (oj/c)  so  that  the  basic  system  (l)-(4)  becomes: 


w^/k  +  uw  +  oau  ^2.  +  uuA j^2A  =  0 

(15) 

u/2  +  cow  /k+  uu  /2  -  9ti  c^/20c  =  0 

t  X  X  'x  V 

(16) 

=  0 

(17) 

where  it  has  been  noted  that 

V®' 


and  equation  (3)  omitted  so  that  the  basic  system  (l)-(4)  is  effectively  reduced 
to  three  equations,  adding  and  substracting  (15)  and  (16)  gives: 

u^/2  +  w^k  +  (u  +  u)[w +u +  uuA j^2A  -  9ri^cV20c^  =  0  (18) 

-  u^/2  +  w^k  +  (u-w)[w 4  ~ 

A  formal  linearization  of  the  equations  (17),  (18)  and  (19)  in  the  neigh¬ 
borhood  of  a  known  isentropic  solution,  denoted  by  the  subscript  zero,  leads 
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to  the  following  system  of  linear  equations  for  the  terms  of  first  order,  denoted 
by  the  subscript  one: 


[w^/k+  Uj/2]^  +  (u^  +  w^)[wQ/k  +  Uq/2]^  +  (u^^  +  a)Q)[Wj/k  +  u^/2]^ 


(20) 


[Wj/k  -  u^/2]^  +  (u^  -  Wj)[wQ/k  -  Uq/2]^  +  (Uq  -  coQ)[w^/k  -  u/2]^ 


(21) 


^It  ^O’^lx  " 


(22) 


According  to  (22),  t|^  remains  constant  along  the  particle  paths  of  the 

given  flow,  i.e.,  along  dx/dt  =  u^  .  Since  PQ(dx  -  u^dt)  is  the  exact 
differential  of  a  function  i}Jq  which  when  equated  to  a  constant,  defines  the 
particle  paths,  the  solution  of  (22)  may  be  written  as: 


*  «(^0^ 

with  12  an  arbitrary  function.  It  is  convenient  to  define  a  new  function 
TQ(x,t)  by: 


(23) 


(24) 
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Introducing  the  functions 

R  =  u^/2  +w^/k  ,  S  =  -  u^/2  +  w^/k 

and  noting  that 

Uq/2  +  WqA  =  q  ,  -  Uq/2  +  w^/k  =  p  , 

the  characteristic  parameters  of  the  basic  flow,  equations  (20)  and  (21) 
may  be  written  as: 

+  (Uq  +  <Jq)R^  +  [3R  -  S  -  (R  +  S)/3Wq^/^]  ay 2 

(25) 

To/2  -  UQtOQA^y2AQ 

St  +  (Uq  -  03^)8^  +  [R  -  3S  +  (R  +  S)/3wQ^/^]py2 

(26) 

=  -  Tq/2  -  UQt0QA^y2AQ 

Equations  (25)  and  (26)  have  proved  useful  in  discussing  perturbations 
of  simple  waves.  For  the  present  problem,  the  base  flow  is  a  constant  state 
so  that  a  and  p  are  constants,  and  the  general  solution  in  terms  of  three 
arbitrary  functions  of  one  argument  is; 

Tq  =  2E'[x  -  Upt]  (27) 

R  »  F[x  -  (Ug  +  WQ)t]  +  E[x  ~  Uptl/w^  -  UQWQA^/2AQ(Ug  +  u^)  (28) 

S  =  G[x  -  (Uq  -  WQ)t]  +  E[x  -  Uptl/cop  -  UqWqA^/2Aq(Uq  -  w^)  (29) 

Thus,  there  are  four  distinct  contributions  to  the  perturbation,  viz., 
a  disturbance  due  directly  to  the  area  variations,  one  due  to  the  entropy 
variations,  which  travels  along  the  particle  paths  and  is  measured  by  E, 
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a  perturbation  propagating  with  velocity  co^  ,  the  true  speed  of  sound,  with 
respect  to  the  fluid  along  the  family  of  characteristics  x  -  (u^  +  =  constant 

and  measured  by  F  and  a  perturbation  propagating  with  velocity  with 

respect  to  the  fluid  along  the  family  of  characteristics  x  -  (u^  -  *  constant 

and  measured  by  G. 


3.  Solution  in  the  vicinity  of  the  incident  shock 

Jump  conditions  across  normal  magnetohydrodynamic  shocks  have  been 
considered  by  several  authors,  e.g.,  Friedrichs  [4],  Let  U  be  the  shock 
velocity,  v  =  U  -  u  and  the  subscripts  one  and  two  designate  flow  quantities 
in  the  regions  in  front  of  and  behind  the  shock.  Then  the  analogs  of  the 
Rankine-Hugoniot  relations  are; 


^I'^l  =  Vz 


PlVi^  +  Pi  +  BiVzix  = 


yP^/(y  -  l)Pj  +  ^1^2  +  BiVp^P'  =  ~  ^Z^^^Z^^ 


For  gas-dynamics  shocks,  knowledge  of  the  flow  in  front  of  and  one  para¬ 
meter  behind  the  shock  suffices  to  give  a  complete  solution.  In  magnetohydro¬ 
dynamics,  all  quantities  behind  the  shock  may  be  expressed  in  terms  of  those 
in  front  and  two  parameters,  viz. ,  the  shock  strength  and  one  which  gives  a 
measure  of  the  applied  field. 

Let  m  =  b/c,  n  =  u/c,  cr  =  P2/P1 »  '''  =  ^z^^l '  ^  ~ 
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2  2  2  7 

q  =  u  /c  =  1  +  m'^  .  Then  the  following  relations  hold,  where  n^  is 
assumed  zero: 


V1/V2  =  P2/P1  =  ==  ~ 


T  =  [4(r  -  1  -  5m ^^(1  -  -  u) 

(31) 

=  [6cr  +  mj^cr(5  +  (r)]/2(4  -  ir) 

(32) 

n2=  (1  -  (r'^)M^(cr/T)^/^ 

=  (o'  -  1)|[3  +  mj^(5  +  (r)/2]/[4(r  -  1 

-  5mj^(l  -  o-)V8]| 

C2Vci^  =  T/,r 

(34) 

m^Vm^^  =  (tVt 

(35) 

2  2  , 

=Mj  /o-T 

(36) 

which  express  the  flow  parameteis  behind  the  shock  in  terms  of  cr  and  m^  . 

The  effect  of  an  area  variation  bn  the  motion  of  an  initially  uniform  shock 
propagating  with  constant  speed  into  a  fluid  at  rest  will  now  be  determined 
by  the  use  of  the  general  solution  for  the  non-isentropic  perturbation  of  a 
constant  state,  equations  (27)-(29).  The  method  of  generation  of  the  shock 

is  left  open  save  that  whatever  that  may  be,  the  mechanism  is  sufficiently 
far  removed  so  that  no  reflections  come  back  to  interfere  with  the  basic 
interaction  considered  herein. 

From  the  formulation  of  the  problem,  there  is  no  mechanism  downstream 
of  the  shock  which  could  give  rise  to  the  term  F[x  -  (u^  -  “2^^^  (28).  The 
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pressure  perturbation  behind  the  shock  will  be  determined  by  setting  F  =  0. 
This  gives: 

u/z  +  w/k  -  .  -  u^co^A/ZA^Cu^  +  co^)  (31) 


where  perturbations  of  a  quantity  are  denoted  by  a  bar.  Since 


w/k  = 


9  n/20c^  =  SP/lOP^  -  3C/2C2  , 


the  latter  from  the  equation  of  state,  equation  (37)  may  be  written  as: 


u/Z  +  /Zq^  +  So^P/lOq^P^  =  -  +  q^jA^ 


Uj/Cj  +  +  ZlPj/lOq^Pj  -  P2/^‘’2P2  °  "  "2'’2*2^*"2  '''  ‘'2*\ 


From  the  jump  conditions,  equations  (30)-(36), 


—  [90(r  +  5m^^cr(l  1  -  24(r  +  15(r^  -  Zcr^)] 

(4  -  cr)[6(4(r  -  1)  -  5m ^^{1  -  o-)^] 


u^  [24  +  IZo"  +  m^  (20  +  lOcr  +  l<y  -  a^)]  n^v 
*^2  2(4  -  cr)(cr  -  1)[6  +  m  ^(5  +  (r)] 


so  that  (38)  gives: 


^2/^2  =  " 
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where: 

r 

n_q_  90or+5m.^o-(l  l-24<r+15(r^-2<r^)  n  24+12<r+m,^(20+10(r+7(r^-(j-^) 

jr  j.  --  ^  ^  -  —  -  _  —  —  ,  ^  ^  _ _ 

^  '^2^*^2  6{4(r-l)-5m^(l-(r)^  2(o--l)[6+mj^(5+o-)] 


3(5m2+2)  90(r+5m^^(r(ll-24(r+15o-^-2o-^) 
^*^^2  6(4o--l)-5m  ^^(1-0-)^ 


(4-0-) 


or 

P2  =  -  mj){P2  -  (40) 

where 

5m^^(l-<r)^  -  6(4(r-l) 

30(l-<r)  +  5tnj^(l-o-)^ 


Strong  shocks  can  occur  in  magnetohydrodynamics  in  two  ways,  viz, , 
for  cr  close  to  4  or  for  a  very  strong  applied  field  for  any  (r>l,  i.  e.,  m^ 
is  then  large. 

For  all  m  ,  Urn  K_  =  0.5  and  lim  K_  *  0.45085,  but  these  are  pre- 
o-*l+  (r-*4- 

cisely  the  limits  for  the  ordinary  gas  dynamic  case  which  corresponds  to  m^  =  0  . 
Thus,  for  very  strong  or  very  weak  shocks  [in  an  asymptotic  sense],  the  results 
are  independent  of  the  magnetic  field  and  agree  with  the  usual  gas  dynamic 
results. 

For  m^  =  0,  K2  decreases  monotonically  with  tr  and  agrees  exactly 
with  the  corresponding  parameter  in  the  gas  dynamic  case  [2], 

For  any  m^  ^  0,  there  no  longer  is  a  monotonic  variation  with  cr,  but 
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the  curves  are  all  concave  upward  with  those  for  greater  lying  beneath 
those  for  lesser  .  All  the  curves  issue  from  the  point  (  o- )  =  (1,0.  5) 
and  terminate  at  the  point  (4,  0.  45085). 

For  fixed  incident  shock  strength  o-,  decreases  monotonically  with 
increasing  m^,  so  that  in  a  diverging  (converging)  channel,  the  pressure 
decrement  (Increment)  is  decreased  (decreased)  by  increasing  the  applied  field. 

Qualitatively,  the  motion  of  the  shock  is  independent  of  whether  the  main 
flow  behind  the  shock  is  <  oo  or  >  w. 

Graphs  of  K2  are  given  in  Figure  1, 


4.  The  Reflected  Disturbance 

To  complete  the  solution,  the  arbitrary  function  G  of  equation  (29)  must 
be  determined.  This  is  done  most  readily  by  noting  that  the  system  (27)- (29) 
may  be  written  as: 

u/2  +  3c2m2V2q2  +  3c2P2/l0^2P2  =  '  ‘  +  ^'^Z^Z^^^Z 

+  ^<^2^Z'^^^^Z^Z  ~  ~  ^z'^Z^Z'^^^Z^^Z  ~  ^2}  ”  ^^2  "  “2^^^ 


Thus,  on  addition: 


15m  ^+6  3m^(4- (t)  [6(4<r-l)-5m  ^(l-cr)^] 


lOq. 


2q  <r 


[90+5m.  (1  l-24(r+15(r  -2cr  )] 

^  A 


Hp  Qp  A 

ir~l  aT  =  ■  '“2  '  “2>'l^‘=2 

"2 ■”2  ^ 
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Evaluating  this  on  the  shock,  x  =  Ut,  and  replacing 
value  from  (39)  gives: 

G[k]/C2  =  K3A2[6k]/A2 


where 


X.  =  X  -  (u^  - 

6  =(M2  * 


2 

2 

2  »3  ,  , 

5m2+2 

m^  (4-0-) 

Sm^  (l-(r)  -  6(4or-l) 

^^2 

5  O' 

90+5mj^(ll-24(r+15(r^-2(r^) 

Thus: 


(41) 


where 


= 


2 


"2“'^2 


5m  j^(l-(r)^-6(4o— 1) 


30(l-v)+5m^^(l-v)^ 


2q-: 


2  2 
5m  +2  m  (4-v) 

- 4 -  +  -± - 


5m^(l-o-)^  -  6(4o— 1) 


I  90+5m^^(ll-24cr+15o-^-2o-^) 


ft 


^2  ■  “  ^1  *^2 

The  parameter  6  is  a  monotonic  increasing  function  of  the  shock  strength 
and  varies  between  the  limits 

0.5  <  6  <  1.  2361 

for  all  m^.  Graphs  are  presented  in  Figure  2. 
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It  might  be  noted  that  Chester's  graph  of  the  corresponding  function  is 
incorrect.  He  gives  it  as  a  monotonic  decreasing  function  and  the  value  for 
a  very  weak  shock  (cr  =  1)  is  the  reciprocal  of  the  true  value. 

When  the  flow  behind  the  shock  is  <  w,  a  disturbance  will  be  reflected 
to  the  left  downstream  of  the  shock,  and,  from  (41),  the  pressure  perturbation 
will  ultimately  be  given  by; 

since  A2  =  0  for  x  <  C  . 

For  flow  >  w  behind  the  incident  shock,  the  reflected  disturbance  will 
travel  to  the  right  and  the  pressure  perturbation  is  given  by  (41). 

Graphs  of  the  parameter  are  presented  in  Figure  3.  The  parameter 
is  singular  for  n^  =  This  difficulty  may  be  eliminated  by  a  more  careful 
consideration  of  the  basic  perturbation  equations  by  retention  of  appropriate 
nonlinear  terms.  However,  since  the  main  purpose  of  this  paper  is  to  consider 
the  pressure  perturbation  in  the  immediate  neighborhood  of  the  incident  shock, 
such  are  not  employed. 


5,  Concluding  Remarks 

In  conclusion,  let  it  be  noted  that  a  theory  has  been  developed  rather 
parallel  to  that  of  Chisnell  [2]  for  the  magnetohydrodynamic  case,  and  this 


will  be  reported  when  the  numerical  computations  have  been  completed. 


For  the  present,  suffice  it  to  say  that  this  theory  gives  that  the  strength 

of  converging  cylindrical  and  spherical  magnetohydrodynamic  shocks  near  the 

-K,  -2K_ 

Zoo  Zoo 

point  of  collapse  is  proportional  to  D  and  D  ,  respectively, 
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independent  of  the  applied  field,  where 
distance  from  the  point  of  collapse. 

The  author  wishes  to  thank  Mr.  J. 
computations  employed  herein. 


K-  =  0.45085  and  D  is  the 
Zoo 


Al-Abdulla  for  carrying  out  the 


Figure  1 


A  graph  of  the  parameter 
vs.  the  shock  strength. 
The  uppermost  curve  corres¬ 
ponds  to  m^  -  0  ,  i.  e. ,  the 
gas  dynamic  case  while  the 
successively  lower  curves 
correspond,  respectively  to 
m ,  =  1 ,  1.  5,  3  and  6. 
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